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Abstract. We present the integrated mass profiles for a sample of ten nearby (z < 0.15), relaxed galaxy clusters, covering a 
temperature range of [2 - 9] keV, observed with XMM-Newton. The mass profiles were derived from the observed gas density 
and temperature profiles under the hypothesis of spherical symmetry and hydrostatic equilibrium. All ten mass profiles are 
well described by an NFW-type profile over the radial range from 0.01 to 0.5i?200> where ^^200 is the radius corresponding 
to a density contrast of 200 with respect to the critical density at the cluster redshift. A King model is inconsistent with 
these data. The derived concentration parameters and total masses are in the range C200 = 4-6 and M200 = 1-2 lO''' - 
1.2 10'^ Mq, respectively. Our qualitative and quantitative study of the mass profile shape shows, for the first time, direct and 
clear observational evidence for the universality of the total mass distribution in clusters. The mass profiles scaled in units 
of R200 and M200 nearly coincide, with a dispersion of less than 15% at 0.1 i?200- The C200-A/200 relation is consistent with 
the predictions of numerical simulations for a ACDM cosmology, taking into account the measurement errors and expected 
intrinsic scatter. Our results provide further strong evidence in favour of the Cold Dark Matter cosmological scenario and show 
that the dark matter collapse is well understood at least down to the cluster scale. 
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1. Introduction 

The Dark Matter (DM) distribution in clusters is a sensitive test 
of current scenarios of structure formation and of the nature 
of the DM itself. Of particular interest are comparisons with 
the predictions of N-body simulations of hierarchical cluster- 
ing in the currently favoured ACDM (Cold Dark Matter) cos- 
mology. Numerical simulations (e.g "Navarro et al. "1997), as 
well as early semi-analytical work (e.g. Bertschinaer 1985), 
predict a remarkable similarity in the Cold Dark Matter density 
profile of virialized halos. Although the exact slope in the very 
centre is still a matter of debate, recent high resolution simula- 
tions predictdia^daAnwtte^ JNavarro et al 
1997* ' Moore et alJll999t iDiema nd et al."2004i: iNavarro et al 



re et al 

2004) and that the concentration of the Dark Matter varies only 
slightly with system mass (e.g.|Dolag^talj2QQi. 

The strong similarity in the ROSAT surface brightness 
profiles iVikhlin ii i et al. 1999; Neumann & Arnaud 2001; 
lArnaud et alJl2002ah . and of the temperatu re profiles of hot 
clusters observed with A5CA and BeppoSAX ('Markevitch et al.' 
Ll998; Irwin & Bregman 2000; De Grandi & Molendi 2002) 
provided indirect evidence of a universal underlying dark mat- 
ter distribution. The present generation of X-ray satellites. 
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XMM-Newton and Chandra, represent a giant step forward in 
terms of resolution and sensitivity. We can now measure pre- 
cisely, through the hydrostatic equilibrium equation, the total 
mass distribution in clusters. Evidence is slowly accumulat- 
ing that CDM numerical simulations predict the correct shape 
of t he Dark Matter distribution, not only in massive clusters 
(e.g.lDavid et alJ200l|lAllen et al.|200lllArabadiis et alJ2002t 
iBuote & Lewisl '2004':'Buote"2004': Allen et al.""2003'). but also 
in low mass clusters cPratt & Ar naud 2003). This may well 
be true up to the viriai radit is, as shown up to 6 = 500 by 
the observation of A1413 jPratt & Arnaudll2002h . The ob- 
served profiles are cusped in the centre, and, for a few mas- 
sive clusters, the inner slope has even been measured pre- 
cisely enou gh to dist inguish be tween various CD M predictions 

jLewis et al 2003: iBuote & Lewis,, .2004: .Pointecouteau et all 

120041) . However, most mass studies ha ve been conducted onin- 
dividual 'test case' clusters. Recently. IPratt & Arnaud! ( 120051) 
performed the first quantitative check of the universality of 
the mass profile using a sample of five clusters observed with 
XMM-Newton (four low mass systems compared to one mas- 
sive system). It is necessary to extend this type of study to 
larger samples, with a better temperature (i.e., mass) coverage. 

In this paper, we use XMM-Newton to examine the total 
mass profile of ten relaxed, nearby clusters in the tempera- 
ture range from 2 to 9 keV. In a companion paper (Arnaud, 



Fointecouteau & Pratt, 2005; Paper 11), we use the mass data 
to study the scaling properties of the mass with temperature. 
In Sect. 121 we present the sample, the observations and the data 
processing steps. We detail the extraction of the scientific prod- 
ucts, from temperature and density profiles to mass profiles. We 
quantify the shape of the mass profiles in Section|3 Our results 
are discussed and we conclude in Sec.|4] 

We have used the currently favoured ACDM cosmology, 
Hq - 70 km/s/Mpc, Q,„ = 0.3 and = 0.7, throughout this 
paper. 

2. XMM-Newton observations and data processing 

2.1. The sample 

The sample has been built to cover a wide range in temperature 
(i.e. in mass) from 2 to 9 keV, and constitutes 10 clusters. We 
limited the redshift range to z < 0.15, where evolution effects 
are expected to be negligible. With the exception of A478 and 
PKS0745, for which there was a mosaic observation, we only 
considered clusters fitting in the XMM-Newton field of view, 
enabling the local background to be estimated, thus limiting 
systematic uncertainties on the temperature, and consequently, 
the mass profiles. Since cluster size increases with tempera- 
ture, this last criterion sets a lower limit on the redshift for 
each temperature. Cool and hot clusters lie in a redshift range 
close to this limit ([0.04,0.06] and [0.1 - 0.15] respectively). 
This thus optimizes both the cluster coverage and the statistical 
quality of the data. A final selection criterion was the quaUty of 
the mass data. All clusters in the sample have a regular X-ray 
morphology, indicative of a relaxed state and allowing reliable 
determination of the total mass profile through the hydrostatic 
equilibrium equation. 

Our sample includes the sample of lPratt & ArnaudI ( l2005l) : 
the cool clusters A1983, A1991, MKW9, A2717, and the 
hot cluster A1413. We improved the temperature coverage by 
adding A478 (recently studied by Fointecouteau et al. 2004 ). 
and four clusters at intermediate and high temperature from the 
XMM-Newton archive which meet our criteria. The journal of 
observations is presented in Tabled 

To minimise systematic errors in the statistical study of 
cluster properties, it is important to use scientific data derived, 
as far as possible, with the same procedure. Unless otherwise 
state d, we use the previously pub lished d ata of the cool clus- 
ters (Pratt & Arnaud'2003', 2005), A478 ( Pointecoutea u et all 
b004) and A1413 ( Pr att & Arnaud 2002)'. These results were 
obtained with the same general method that we use to process 
the four additional clusters, although some details differ from 
cluster to cluster. The procedure is described in the next four 
sections. Further details and comments on each individual tar- 
get are given in AppendixlAl 

2.2. Event list processing 

We made use of the XMM-Newton sas software package, 
versions 5.3 or 5.4, to filter the data. Below we detail the main 

' We rescaled the published values of A 1983 and A 14 13 to the 
ACDM cosmology used in this paper. 
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Notes: " EPN observation mode: FF= Full Frame, EFF=Extended Full 
Frame; ^ Exposure time (EMOS1/EMOS2/EPN) in ksec after flare 
cleaning; The exposure times are given for the central pointing. 

data processing steps. 

1 . Considering only events with FLAG - and PATTERN < 
12 (EMOS) and PATTERN = (EPN), we clean the 
data for soft protor i flares using a thresh old cut method 
(see appendix A in lPratt & Arnaudll2002l) . A first screen- 
ing was performed using light curves in the high energy 
band ([10-12] keV for EMOS, [12-14] keV for EPN) us- 
ing 100s bins. After filtering using the good time intervals 
from this screening, a second screening was performed as a 
safety check-up. A second light curve built in 10s bins in a 
wider energy band ([0.3-10] keV for EMOS, [0.3-12] keV 
for EPN) was screened, and the event lists were filtered ac- 
cordingly. 

2. To coiTect for the vignetting effect, w e used the photon 
weighting method jArnaud et al.l200ll) . The weight coeffi- 
cients were computed by applying the sas task evigweight 
to each event file (including background files). 

3. The point source lists from the SOC pipeline were visually 
checked on images generated for each detector Selected 
point sources from all available detectors were gathered 
into a single point source list and the events in the corre- 
sponding regions were removed from the event lists. The 
area lost due to point source exclusion (as well as CCD 
gaps and bad pixels) was computed using a mask image. 

We used XMM-Newton de dicated blank field datasets from ei- 
flier lLumb etaP ( l2002l) . or lRead & PonmanI ( l2003l) to obtain a 
background event list associated with each data set^. For each 
cluster, the blank field event list was recast in order to match 
the astrometry of the observation and was then processed in 
the same way as the observation event list. (This included ex- 
traction of the events in the same regions as the point sources 
in the observation data set.) To account for variations in the 
particle background level, we assumed that only particle sig- 
nal is collected at high energy, and used the ratio between the 

^ Most of observations used for the lRead & Ponmar] i2003h blank 
field were centred on point sources. Their removal, and the removal of 
other point sources in each field, produced local variations of a factor 
of two from point to point in the final stacked exposure map. We took 
this variation into account in our analysis. 



high energy count rate or the observation and the blank held to 
renormalise the blank field counts. The bands used were [10- 
12] keV for EMOS and [12-14] keV for EPN. To avoid the 
contamination due to possible high energy emission from the 
cluster, a central region (r < 5') was excluded when comput- 
ing the normalisation in the case of the hot clusters. For EPN, 
the out-of-time (OoT) events were considered as an additional 
background component. For each observation the OoT event 
list was generated using the sas task epchain and processed in 
the same way as the observation. A normalisation factor was 
then used depending on the observing mode (Full Frame or 
Extended Full Frame). 

The cleaned event list, the blank field and, in the case of the 
EPN, the OoT event file of each single observation were used 
to extract scientific products such as spectra and surface bright- 
ness profiles. Background subtraction was per formed using 
the dou ble subtraction process fully described in lArnaud et al.l 
Appendix), and involves subtraction of the normalised 
blank field data, and subsequent subtraction of the Cosmic X- 
ray background residual estimated from a region free of cluster 
emission. 

2.3. The density profile 

For each cluster, an azimuthally-averaged, background sub- 
tracted surface brightness (SB) profile was computed in the soft 
energy band ([0.3-2.] keV in the present work) for each avail- 
able detector^. The profiles of all detectors were then summed 
together into a total SB profile and rebinned with logarithmic 
radial binning and a minimum S/N ratio of 3cr. The profile 
was corrected for radial variations of the emissivity (e.g due to 
abundance or temperature gradients) in the energy band con- 
sidered (see Pratt & Arnaud 2003, for details). This corrected 
profile is thus proportional to the emission measure along the 
line of sight. 

The final SB profiles were fitted using parametric ana- 
lytic models of the gas density profile, converted to an emis- 
sion measure profile and convolved with the PSF spatial re- 
sponse CGhizzard i 2001, 2002). We considered various para- 
metric forms and empirically chose the model best fitting the 
data using the statistic as a measure of the goodness of fit. 
T he models included a d ouble y6-model (the BB model defined 
m IPratt & Arna'u3i2002l) . a modified double j8-model, which 
allows a more concentrated gas density distribution towards 
the centre (the KBB model used by Pratt & Arnaud 2002 for 
A1413), and the sum of three jS-models in which a common 
value of /3 is assumed to ensure smooth behaviour at large radii 
(the BBB model used by Pointecouteau et al. 2004 to model 
the SB profile of A478). For the clusters newly analysed here, 
the best fitting model was either a KBB model (A2597, A1068 
and PKS 0745-191) or a BBB model (A2204). Over the whole 
sample the reduced values vary from 1.1 (for A478), to 1.5 

We also computed background subtracted SB profiles in the high 
energy band (we recall: [10-12] keV for EMOS and [12-14] keV for 
EPN) to check that there were no significant residuals, i.e., that the 
normalisation of the blank field background was correct, and that it 
was indeed constant with radius. 



(tor AlOos), and thus even the best parametric model leads in 
some cases to a formally unacceptable fit. This is linked to the 
very small statistical errors on each measurement. The actual 
discrepancies between the model and the data remain small: 
adding a few percent (2 to 5%) of systematic error while fitting 
the data always leads to acceptable values. 

2.4. The temperature profile 

Concentric annular regions were defined from the total back- 
ground subtracted SB profile of the cluster the clusters newly 
analysed here, we used the following empirical method to de- 
fine the annuli from which the spectra were extracted. We 
started with a minimum bin size of 15" (e.g about the HEW 
of the XMM-Newton PSF) and then increase the bin size by a 
logarithmic factor of 1.05. This is sufficient in the central parts 
to keep a good S/N, but not in the outer regions. We thus further 
impose, in the outer regions, that the numbers of cluster counts 
per bin is approximately constant, within 5%. 

For each annular region, a background subtracted spec- 
trum was extracted for each available detector The spectra 
from e ach detector were then simultaneously fitted using xspec 
(lArnaud.1996.) with an absorbed, redshifted single tempera- 
ture plasmaniodelJ^BS^[\ffiKAL). Except in the case of A478 
(see lPointecouteau et alJ200 4][ for details), after checking that 
the A^H value agreed with the galactic 21 cm value (from 
iHartm ann & Burton lll999l) '). this parameter was frozen. Thus 
we derived an annular temperature profile for each cluster in 
the sample, as shown in the left panel of Fig. ^ All profiles 
(except that of A1413) have the same generic shape: a negative 
gradient towards the centre and a roughly flat external plateau. 
A1413 is the only cluster for which we observe a significantly 
decreasing temperature at high radii (of about 20%). 

To derive the true radial temperature profile (needed to 
compute the mass profile) we should take into account projec- 
tion and PSF effects. For the typical temperature profile shape 
we have obtained, these effects depend most strongly on the 
magnitude of the gradient in the centre. In addition, a strong 
temperature gradient in the central regions is usually associ- 
ated with a strongly peaked surface brightness profile, further 
increasing the PSF blurring^. Reconstructing the radial temper- 
ature profile in these cases is not a trivial task. General non- 
parametric methods, such as simultaneous fitting of annular 
spectra, amplify the noise considerably. This yields radial tem- 
perature profiles with unphysically large fluctuations, particu- 
larly when bot h PSF and proje ction effects are important (see 
discussion in P ointecouteau et al.-2004) . 

The clusters newly analysed in the present work have quite 
strong 'cooli ng flows'. We thus applied t he method developed 
for A478 by IPointecouteau et alJ (l2004 . in which the radial 
temperature profile is derived from the annular temperature 
profile in the following manner. The noise amplification prob- 
lem is avoided by using smooth parametric representation of 
the annular temperature profile. We used the function given by 
lAllen et al.l 1I2OOII) : T(r) ^ To + Ti [ ^^^^ \ to fit the an- 

^ The effects of the PSF correction are most noticeable for A2204 
and A478. 
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Fig. 1. Left: Annular (projected) cluster temperature profiles. Right: Deprojected, PSF-corrected temperature profiles. 
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nular profile. The best fitting model profile is then corrected 
for both the projection and PSF effects, assuming that the an- 
nular temperatures are em ission weighted temperatures (see 
[Pointecouteau et af] l2004 for details). To estimate the errors 
we repeated the procedure 1000 times, using a Monte Carlo 
method that randomizes the annular profile based on the ob- 
served errors. Because using a specific functional form effec- 
tively limits the allowed profiles, the standard deviation of the 
corrected temperature at a given point is occasionally smaller 
than the error on the annular temperature. When this was the 
case, we kept the observed error. 

For the poor systems, which have a more modest cen- 
tral temperature gradient, the radial temp e rature profile was 
estimated as described in IPratt & ArnaudI ( l2003l 12005). For 
A1413 PSF and projection effects proved to be negligible 
jPratt & Arnaiidl E002). The deprojected, PSF corrected pro- 
files of all the clusters are shown in the right panel of Fig.[n 
A detailed discussion of the shape of these temperature profiles 
is beyond the scope of this paper. 



2.5. The total mass profile 

The mass profile for each cluster was derived from the best fit- 
ting density profile and the deprojected, PSF corrected temper- 
ature profile under the assumptions of hydrostatic equilibrium 
and spherical symmetry: 



M(r) = - 



kT r 
Gfinip 



dlnria dlnT 
+ 



dlnr dlnr 



(1) 



The mass, and associated errors, were then computed at each 
radius of the temperature profile^ using a Monte Carlo method 
that ran domises the temp erature profile based on the observed 
errors jPrattfe ArnaudI I2OO3I) . A cubic spline interpolation 

^ We use in fa ct the weighted effective radius of each annulus 
jL^isitai]2QQ3). 



is used to compute the temperature derivative at each radius. 
Only mass profiles which monotonically increase with radius 
were kept (1000 in total). The errors due to the errors on the 
density gradients, which are negligible with respect to those 
from the temperature, were then added. The mass profile of 
A 14 13 was originally derived using a slightly different Monte 
Carlo method (Pratt & Arnaud 2002). For homogeneity of the 
present study, we re-derived the mass profile of this cluster with 
the present method, allowing us to extend the profile to lower 
radii^. The new mass profile is perfectly consistent with the 
published one. 

Each cluster has been checked for the presence of struc- 
ture such as cold fronts, hot bubbles, ghost cavities, or other 
effects which could affect the mass determination. Details of 
each cluster are given in AppendixlAl 

3. The shape of the mass profile 

3.1. Mass profile modelling 

For each cluster, three mass models were fitted to the data: 
(i) a King isothermal sphe re profile; (ii) a standard NFW 
profile dN avarro e t all 1 1 9971) : and (iii) an MQGSL profile 
d Moore et: aljl999i) . Our data indicate that an isothermal sphere 
model (i.e., a profile with a core) is not a good representation of 
the mass distribution in these clusters. Dropping too rapidly in 
the centre and flattening in the outer regions, the reduced ob- 
tained from King model fits ranged from 1.65 for A478 to 15.6 



* IPratt & Arnau j i2002|l use d the original Monte Carlo method of 
iNeumann & Bohrinserl ilQQsl) . This method uses a 'diffusive' pro- 
cess to calculate random temperatures at equally spaced radii. It has 
two parameters, the radial step and a window parameter, which are op- 
timised to ensure smooth random temperature profiles. In this original 
version of the code, the radial step fixed the minimum radius of the 
mass profile. In the present method it is simply the radius of the first 
temperature annulus. 
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Fig. 2. Left panel: Integrated total mass profiles plotted in units of physical radius (kpc). The solid lines are the best fitting NFW profiles as 
detailed in Table|2| Right panel: Scaled mass profiles of all clusters. The mass is scaled to M200, and the radius to R200, both values being derived 
from the best fitting NFW model. The solid black line corresponds to the mean scaled NFW profile and the two dashed lines are the associated 
standard deviation. 



Table 2. Results for the NFW mass profile fits. 



Cluster 


C200 


^^200 (kpc) 


MzoodO'^Me) 


R500 (kpc) 


MjoodO'^Mo) 


^^/^wCdof) 




A1983 


3.83 ±0.71 


1 100 ± 140 


1.59 ±0.61 


717 ±79 


1.09 ± 0.37 


9.1(7) 


[31, 12] 


MKW9 


5.41 ± 0.67 


1006 ± 84 


1.20 ±0.30 


668 ±51 


0.88 ± 0.20 


3.6(8) 


[24, 5.2] 


A2717 


4.21 ± 0.25 


1096 ± 44 


1.57 ±0.19 


717 ±26 


1.10±0.12 


16(10) 


[48, 52] 


A1991 


5.78 ±0.35 


1106 ±41 


1.63 ±0.18 


737 ± 25 


1.20 ±0.12 


10(9) 


[82, 33] 


A2597 


5.86 ±0.50 


1344 ±49 


3.00 ±0.33 


897 ± 29 


2.22 ± 0.22 


14.6(8) 


[55, 6.8] 


A1068 


3.69 ± 0.26 


1635 ±47 


5.68 ± 0.49 


1060 ± 26 


3.87 ± 0.28 


2.5(4) 


[62, 9.8] 


A1413 


5.82 ±0.50 


1707 ± 57 


6.50 ± 0.65 


1129 ±33 


4.82 ± 0.42 


8.0(8) 


[45, 3.7] 


A478 


4.22 ± 0.39 


2060 ± 110 


10.8 ± 1.8 


1348 ± 64 


7.57 ± 1.11 


9.5(10) 


[16, 30] 


PKS0745 


5.12 ±0.40 


1999 ± 77 


10.0 ± 1.2 


1323 ± 45 


7.27 ± 0.75 


2.3(6) 


[20, 20] 


A2204 


4.59 ± 0.37 


2075 ± 77 


11.8± 1.3 


1365 ± 44 


8.39 ±0.81 


9.7(6) 


[14, 21] 



(a) Chi-square obtained for the best fit of a king and a MQGSL model respectively (see Sect. l3.1> 



for A 1068. It is rejected with a minimum 91% confidence level 
(A478). In contrast, the reduced obtained from NFW model 
fits varied from from 0.4 (PKS 0745-191)to 1.8 (A2597), while 
the MQGSL profile yielded reduced of 0.5 (A1413) to 5.21 
(A2717). We note that the chi-squared value is very sensitive 
to the central points. In some cases, the mass errors on these 
points may be underestimated due to the procedure used for 
PSF and projection effects correction, or to systematic errors 
we are not able to quantify and are therefore unable to take into 
account. 

In all cases barring A1413 jPratt & ArnaudI l2002h and 
A2597, the NFW profile proved to be a better fit than the 
MQGSL model, and de facto to be the best representation of 
our current data. In the cases of A1413 and A2597, the im- 
provement inx^ when an MQGSL profile is used is very small. 
Thus, to keep our approach coherent, we decided to use the 
NFW fit as a parametric representation of the mass profile of 
each cluster. 



The NFW model, where the density is p(r) oc [(r/rs)(l + 
'"/'"s)] has two free parameters: the scaling radius rs, and a 
normalisation parameter. The model can be equivalently ex- 
pressed in terms of the concentration parameter C200 - '"20{)/'"s 
and the total mass Mioo. M200 is the mass corresponding to 
a density contrast of 6 = 200, i.e. the mass contained in a 
sphere of radius R200, which encompasses a mean density of 
200 times the critical density at the cluster redshift: pdz) - 
3E(zfHo^ /SttG, where E^(z) = Qm(l + zf + Oa. In numerical 
simulations, this sphere is found to correspond roughly to the 
virialised part of clusters. The results of the best NFW fits are 
detailed in Table |2l and the best fitting profiles are shown in 
Fig. El 

3.2. Scaled mass profiles 

The left panel of Fig.|2]shows the mass profiles (and best-fitting 
NFW models) plotted in physical units (kpc). Not surprisingly, 
there is a continuous increase in mass with cluster tempera- 




2 - 



the concentration parameter and the cluster mass usmg our ex- 
te nded sample. Figure |3]presents the C200 - A^200 relation found 
bv lDola^ etai] ll2004 (1 +z)c = co(M/Mo)" with a ~ -0.1 in a 
ACDM cosmology with erg = 0.9. All ten points agree with the 
predicted relation, taking into account its intrinsic dispersion 
and the measurement errors. 

We then performed a linear regression fit in the log C200- 
logMaoo plane, taking into account the uncertainties on each 
quantity. The resulting slope of a - -0.04 + 0.03 (x^ (dof) - 
42(8)) is poorly constrained, due in part to the limited size of 
the sample and the large relative uncertainties in the concen- 
tration parameter. The result is however compatible with the 
intrinsic dispersion of theoretical predictions. The best fitting 
c(M) relation is shown as the long dashed hne in Fig.|3] 



Fig. 3. Concentration parameter C200 versus the cluster mass M^oo- 
The solid line represents the variation of C200 for clusters at z = from 
the numerical simulations of lPolas et alJ i2004) . The dotted lines are 
the standard deviation associated with this relation. The dashed line 
represents the same relation at a redshift of z = 0.1 5( the maximum 
redshift for our sample). The long-dashed line stands for our best fit 
of the C200 - ^200 (see text). 



ture, reflecting the temperature coverage of the sample. These 
unsealed mass profiles already show signs of an underlying 
similarity in the matter distribution. The right panel of Fig. |2 
shows the scaled mass profiles, where we express the radius in 
terms of 7?2oo and the mass in terms of M200, these values being 
derived from the best fitting NFW model of each cluster. The 
scaled mass profiles cover a wide range of radii, from about 
0.01 /?200 to 0.7 R200, and are particularly well constrained be- 
tween 0.1 /?200 and 0.5 Rioo- The agreement between the scaled 
profiles is remarkable, reflecting the similarity in shape of the 
profiles. The average of all best fitting NFW models is shown 
as a black line, with dashed lines representing the mean plus 
or minus the standard deviation. The dispersion is small and 
virtually identical to that derived bv IPratt & Arnaud ( l2005l) 
from their smaller sample: we obtained a dispersion of 7.4% 
at 0.3/^200 and 14.3% at O.lT^aw (compared to the 8% and 15% 
found by those authors). 

3.3. Variation of tine concentration parameter witli 
mass 

Structure formation models do not in fact predict a strictly 
universal matter distribution in clusters. A weak variation in 
the concentration is expected from low to high mass clus- 
ters, reflecting differences in the formation epochs of low and 
high mass haloes (Navarro et al. 1997; Bullock etal. 2001; 
Ibolag etani2004l) . Building on the work of IPratt & ArnaudI 
ll2005[ their Fig. 12), we can investigate the relation between 



4. Discussion and conclusion 

We have measured the total mass profile of ten clusters from 
0.01/?200 up to 0.5/? 200 using XMM-Newton observations. Our 
sample has an excellent temperature coverage and covers an 
order of magnitude in mass from M2()o = 1.2 10"* Mq to 
1.2 10'^ Mq. Our study confirms previous XMM-Newton and 
Chandra studies conducted on individual targets (see Sec. 
and extends the initial statis t ical st udy of cluster mass profile 
structure bv lPratt & Arnau 3 il200.5ll . 

We have found that the NFW profile is a good representa- 
tion of the ten observed mass profiles, and that in all cases the 
isothermal sphere model (i.e a profile with a core) is rejected at 
high confidence. In other words, we confirm the cusped nature 
of the Dark Matter profile, as predicted by CDM simulations 
of hierarchical structure formation, over the temperature/mass 
range of the present sample. The mass profile shape is close 
to universal, again as predicted, with a dispersion of less than 
15% at O.lRova in the scaled mass profiles. The shape is quan- 
titatively consistent with theoretical predictions. The variation 
of the observed concentration parameters with mass is in line 
with the predictions, taking into account the measurement er- 
rors and the expected intrinsic scatter. However, our sample is 
still too small to draw any firm conclusions on the exact form 
of the c(M) relation. Taken together, our results provide fur- 
ther strong evidence in favour of the Cold Dark Matter cosmo- 
logical scenario, and show that the physics of the Dark Matter 
collapse is well understood. 

The exact inner slope of the CDM distribution in haloes 
remains an importan t theoretical issue (biemand et alJl2004t 
iNavarro et al.ll2004l) . Few of our mass profiles have the re- 
quired radial coverage and statistical quality in the central parts 
to allow us to firmly distinguish between an NFW-type profile 
and other types of cusped DM profile (e.g., the mass profile of 
A478). We caution also that the very central parts of clusters 
are regions of complex phenomena (hot bubbles, ghost cavi- 
ties, cold fronts, interaction with the central galaxy, etc) which 
are still not well understood. Their effect on the ICM may chal- 
lenge the hypothesis of hydrostatic equilibrium, and therefore 
the reliability of the X-ray mass estimate. 

Observations are also needed of the outskirts of clusters. 
To date, those regions are basically unknown to observers, and 
can only be investigated with numerical simulations. Study of 



these regions is clearly needed to advance our understanding or 
structure formation and evolution. 
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Appendix A: Cluster specifics 

For each cluster in our sample, we searched the literature for 
special features such as ICM bubbles, ghost cavities, cold 
fronts or other phenomena which could disturb the relaxed 
structure of the intra-cluster medium and therefore invali- 
date out hypothesis of spherical symmetry and hydrostatic 
equilibrium. We particularly investigated existing Chandra 
observations of central substructure in clusters, capitalising 
on its high spatial resolution. We used this information to 
check our mass profiles and if needed to exclude one (or more) 
point(s) containing such substructure. Details and notes are 
now given for individual clusters. 

A478 - From the Chandra observations. ISun et al.l ll2003h re- 
ported the presence of X-ray cavities within the 15 cen- 
tral kpc. The recent XMM -Newton observation of A478 
dPointecouteau et ani2Q04l) excluded this inner area for the 
mass profile computation. We used the pubUshed mass pro- 
file of Pointecouteau et al. ( 2004). 

A 1068 - We chose this cluster for its intermediate tempera- 
ture of about 4-5 keV. The data for A 1068 were retrieved 
from the XMM-Newto n archive. A c cordin g to the Chandra 
observation detailed in lWise et alJ li2004l) . no special inner 
structure has been seen. 

A 1 4 13 - This cl u ster h as been studied in detail by 
IPratt & ArnaudI ( l2002h . Among the sample it is the only 
hot cluster without a cool core, and the only cluster with a 
firmly observed decrease in temperature at large radii. It is 
also the cluster observed up to the largest radius (0.7/?20{))- 
The existing (public) Chandra observation of this cluster 
has a too low an exposure time to show clearly whether 
there is structure in the cluster centre. 

A1983 - IPratt & ArnaudI ll2003l) published an extensive study 
of the XMM-Newton observations for this cluster We used 
all the mass points published by those authors. No Chandra 
X-ray observation is yet reported for this cluster. 

A1991 - Sharma et al. (2004), using Chandra, have reported 
low energy knots (kT ~ 0.8 keV) located at ~ 10 kpc 
from the cluster centre. The first bin of the temperature 
(and thus mass) profile, used by Pratt & Arnaud (2005) in 
their XMM-Newton study of A1991, extends up to 38 arc- 
sec. The knot structure seen by Chandra is included in this 
bin and is therefore diluted within the cluster emission. 
Keeping or removing the first point in our mass analysis 
leads to mass modelling results which are compatible at a 



Icr level. Ihe NFW fit being slightly better when using all 
the point s, we choose to use the whole mass profile pub- 
lished bv lPratt & ArnaudI ll2005l) . 

A2204 - This is a very well relaxed cluster at a redshift of 
z = 0.1523 showing a very peaked emission profile toward 
the centre. The cluster hosts a very s t rong r adio source at 
its centre, and recentlv fSanders et al.l (|2004') have pointed 
out complex structure in the very central part of the clus- 
ter using the Chandra satellite. In addition, as discovered 
in many other clusters, two cold fronts were observed at 
radii of about 30 and 80 kpc. While we computed the mass 
profile down to small radii for this cluster, we excluded in 
the mass analysis the two inner annuli, the outer radius of 
which extend up to 109 kpc, in order to avoid invalidating 
the assumption of hydrostatic equ ihbrium. 

A2717 - Another poor cluster from'P ratt & ArnaudI ll2005l) . 
No high resolution X-ray imaging has been performed 
to date on this cluster. However, a work based on 
ROSAT/PSPC and APM optical data aiming to find sub- 
structure in a sample of clusters (Kolokotronis et al. 20Q l|). 
reported substructure in the optical for A2717. It is unde- 
tectable in this X-ray observation. No significant change 
occurring when excluding the central point from the fit, we 
used the entire computed mass profile. 

A2597 - This cluster was taken from the XMM-Newton 
archive, and was chosen because of its intermediate tem- 
perature of about 3-4 keV. A2597 has been studied with the 
Chandra satellite, exhibiting ghost cavities in its centre at a 
radius of about 30 kpc (McNamara et al. 2001). Moreover, 
the Chandra SB profile and temperature profile seems to be 
quite disturbed up to a few tens of kpc. We therefore choose 
to ignore the two first data points of our mass profile and to 
only keep points with inner radii larger than 50". 

MKW9 - This cluster is the least regular of the three cool sys- 
tems recently studied by Pratt & Arnaud (.2005.) . However, 
lacking any external information on its central structure we 
use their published mass profile. 

PKS 0745-191 - A previous analysis of XMM-Newton data 
was published bv lChen et al.l (l2003l) . However, we decided 
to reprocess the whole dataset to keep the data processing in 
our sample coherent. The data for this cluster are strongly 
disturbed by important flaring periods. After flare cleaning, 
we kept only the MOS data, the EPN data being still very 
noisy. Our temper a ture p rofile is in total agreement with 
that of IChen et alJ (l2003l) . a nd also with the measuremen t 
by BeppoSAX at large radii (' Pe Grandi & Molendil 19991) . 
Note that the overall normalisation of our temperature 
profile is somehow lower to that obtained by iHicks et all 
( I2OO2I) with the Chandra satellite. Indeed, these authors 
obtained a plateau temperature of about 10.5 keV above 
200 kpc from the centre against our value of ~ 8 keV. 
Regarding the SB profile from the Chandra observation, 
which exhibits a sharp central cusp (mainly due to the cD 
galaxy), we improved the fit of the mass profile by exclud- 
ing our central bin which extends up to 34". 
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